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A  Model  of  2  Dimensional  Flow  in  a 
Twin  Screw  Mix-Extruder  Kneading  Disc 
Using  a  Full  2  Dimensional  Solution  of 
the  Flow  Equations 


1  Introduction 

The  study  of  the  flow  behaviour  of  propellant  materials  in  twin  screw  mix- 
extruders  (SMEsl  continues  to  be  a  subject  of  considerable  interest  in  the 
international  propellant  community.  Modelling  of  the  flow  with  a  computer 
offers  the  possibility  of  estimating  the  effect  of  various  processing  parameters 
without  the  need  for  many  expensive  experimental  trials.  Recent  work  in 
computer  modelling  of  flow  in  SMEs  has  been  reviewed  in  reference  1,  which  also 
presented  a  simplified  model  of  flow  in  a  2  dimensional  kneading  disc.  The  basis 
of  the  model  was  the  FAN  approach  of  Tadmor  [2],  which  essentially  reduces  the 
dimensionality  of  the  model  by  1,  so  that  the  flow  in  the  2  dimensional  kneading 
disc  was  modelled  by  1  dimensional  flow,  and  flow  in  a  3  dimensional  screw 
channel  would  be  modelled  as  2  dimensional  flow.  While  the  FAN  approach  has 
the  advantage  of  reducing  the  size  of  the  computational  problem,  it  ignores  the 
circulating  component  of  the  flow  which  is  important  in  determining  degree  of 
mixing.  Hence  for  the  study  of  mixing  in  3  dimensional  screw  channels  a  fully  3 
dimensional  model  is  required. 

This  paper  presents  a  model  of  flow  which  does  not  involve  an  artificial 
reduction  of  the  dimensionality  of  the  flow.  The  model  is  evaluated  for  the  case  of 
full  2  dimensional  flow  in  a  2  dimensional  kneading  disc  as  a  first  step,  and  in 
future  it  is  expected  that  the  model  will  be  expanded  to  handle  3  dimensional 
problems. 


2.  Theory 


2.1  Derivation  of  Flow  Equations 

The  flow  in  the  2  dimensional  kneading  disc  is  illustrated  in  Figure  1.  As  was 
done  previously  with  the  FAN  based  model  [1],  the  flow  channel  between  the 
extruder  wall  and  the  kneading  disc  has  been  flattened  out,  with  the  barrel  wall 
becoming  a  straight  line  moving  to  the  right  with  velocity  V.  Since  in  future  the 
model  is  to  be  extended  to  3  dimensions,  where  the  1  direction  will  be  along  the 
flow  channel,  the  cross  channel  directions  considered  here  will  be  taken  to  be  the  2 
and  3  directions.  In  the  3  direction  the  flow  is  solely  pressure  generated,  as  there 
are  no  external  applied  velocities  at  the  boundaries  in  the  3  direction  .  However, 
in  the  2  direction  a  velocity  is  applied  to  the  top  surface  of  the  fluid,  and  hence  the 
velocity  in  the  2  direction  of  the  fluid  in  the  interior  of  the  channel  is,  in  general,  a 
complex  function  of  drag  and  pressure  flow.  However,  the  flow  will  be  assumed 
to  be  locally  Newtonian,  and  therefore  linear,  so  the  effect  of  the  applied  velocity 
can  be  simplified  by  considering  the  flow  to  be  a  linear  combination  of  pressure 
driven  flow  and  imposed  shear  flow. 


0  2 


Kneading  Disc 


Figure  1:  Geometry  of  the  flow  channel  and  mesh  of  cells  used  in  solution  of  flow 
equations. 


The  basis  of  the  model  is  to  divide  the  flow  region  into  rectangular  cells.  Within 
each  cell  the  flow  will  be  assumed  to  be  locally  Newtonian  with  a  constant 
viscosity  in  the  cell,  but  a  different  viscosity  in  each  cell  which  depends  on  the 
average  shear  rate  in  the  cell.  The  rheological  behaviour  of  the  fluid  in  the  flow 
channel  will  be  modelled  using  a  generalised  Newtonian  equation:- 


T  =  HY 


Where  t  is  the  shear  stress,  |i  is  the  generalised  viscosity,  and  7  is  the  shear  rate. 
The  viscosity  will  be  given  by  a  power  law  equation:- 

(i  =  k.  y"’1 

The  equations  of  motion  for  2  dimensional  creeping  flow  of  an  incompressible 
fluid  in  the  2  and  3  directions  are:- 
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These  are  coupled  partial  differential  equations  which  cannot  be  solved  directly. 

The  equations  are  solved  by  producing  an  estimate  of  the  pressure,  and  then 
using  the  pressure  distribution  to  calculate  the  velocity  distribution.  The 
velocities  so  obtained  are  used  to  produce  a  new  estimate  of  pressure,  and  the 
process  is  repeated  until  convergence  is  obtained.  The  details  of  the  solution 
methods  will  be  different  for  pressure  and  velocity. 

The  method  for  solving  for  pressure  is  based  on  a  generalisation  of  the  FAN 
method  used  previously  [2].  The  flow  in  each  cell  can  be  considered  to  be  made 
up  of  pressure  driven  flow  in  the  interior  of  the  cell  with  superposed  shear  flow 
from  moving  boundaries,  as  illustrated  for  flow  in  the  2  direction  in  Figure  2.  The 
values  of  v2t  and  v2b  are  calculated  from  the  total  flow  in  the  channel.  It  is 
important  to  note  that  in  the  model  the  pressure  will  be  evaluated  at  the  centre  of 
the  cell. 

Within  the  cell  the  viscosity  is  assumed  to  be  constant,  and  so  in  evaluating  the 
double  derivatives 
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the  viscosity  can  be  taken  outside  the  derivative.  However,  the  derivatives 
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will  be  evaluated  between  adjacent  cells,  and  so  the  viscosity  has  to  be  treated  as  a 
variable  when  expanding  these  derivatives. 
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Figure  2:  Schematic  velocity  profile  in  2  direction  at  the  boundary  between  cells  (J,  K) 
and  1+1,  K). 


Using  these  assumptions,  equations  1  and  2  can  be  rearranged  to  give:- 
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A  simple  way  to  solve  these  equations  is  by  an  iterative  procedure  where  the 
velocity  field  from  the  previous  iteration  is  used  to  calculate  the  derivatives:- 
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These  derivatives  can  then  be  treated  as  constants  and  be  lumped  with  the 
pressure  derivatives,  and  the  equations  integrated  directly. 


For  the  geometry  illustrated  in  Figure  1,  the  equations  can  be  integrated  to  give 
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where  DX2  and  DX3  are  the  cell  dimensions  in  the  2  and  3  directions. 

The  flows  through  the  boundaries  of  the  cell  in  each  direction  are  obtained  by 
integrating  equations  5  and  6  to  give:- 
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If  the  flow  in  the  positive  2  direction  out  of  cell  (J,K)  is  given  by  Q2(J,K),  and  the 
flow  into  the  cell  by  Q2(J-1,K),  and  in  the  3  direction  by  Q3(J/K)  and  Q3(J,K-1), 
then  conservation  of  mass  requires  that 

Q2(J,K)  +  Q3(J,K)  -  Q2(J-1,K)  -  Q30,K-1)  =  0  (9) 

If  equailuns  7  and  8  for  the  appropriate  cc'1®  are  substituted  into  9,  the  pressure  in 
the  cell  is  given  in  terms  of  the  flow  velocities  and  viscosities.  If  values  of  the 
velocities  and  viscosities  are  available  from  a  previous  iteration,  they  can  be 
substituted  into  conservation  equation  to  give  an  equation  for  the  pressure  in  each 
cell  in  terms  of  the  pressures  in  the  surrounding  cells.  The  method  of  solving  this 
set  of  equations  for  the  pressures  in  each  of  the  cells  will  be  discussed  in  the  next 
section. 

Once  an  estimate  of  the  pressure  distribution  is  known,  it  can  be  used  to 
calculate  an  estimate  of  the  velocity  distribution  from  equations  3  and  4. 

The  method  of  solution  of  these  equations  will  be  discussed  in  a  subsequent 
section.  This  procedure  can  be  iterated  until  the  values  of  pressure  and  velocity 
converge  to  constant  values. 


2.2  Method  of  Solution  of  Equations  for  Pressure  Distribution 


The  pressure  distribution  is  calculated  by  converting  the  differential  equations  7 
and  8  into  finite  difference  equations,  and  substituting  them  into  equation  9.  The 
finite  difference  equation  for  Q2  is:- 
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where  V2(J,K)  is  the  velocity  in  the  2  direction  at  the  lower  right  comer  of  cell  Q,K) 
and  the  subscripted  viscosity  terms  are  the  averages  of  viscosity  at  the  centre 
points  of  the  walls  dividing  the  cells  in  the  2  and  3  directions  respectively.  The 
second  derivative  in  the  2  direction  can  only  be  evaluated  at  the  cell  walls  because 
that  is  where  the  velocities  are  defined.  The  required  value  is  at  the  cell  centre,  so 
the  average  of  the  wall  values  is  used.  If  equation  10  is  substituted  into  equation  9 
and  rearranged  to  put  the  pressure  terms  on  the  left  hand  side  and  other  terms  on 
the  right  hand  side,  then  the  contribution  from  Q2(J,K)  is  given  by:- 
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A  similar  equation  is  obtained  for  Q3(J,K)> 
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The  pressure  term  P(J,K  +  1)  has  been  put  on  the  right  hand  side  so  that  the  left 
hand  side  only  contains  variable  terms  in  K,  and  not  K-l  or  K+l. 

Similar  equations  are  obtained  for  the  flows  Q2(J-1,K)  and  Q3(J,K-1 ). 
Substituting  into  equation  13  gives  an  equation  of  the  form 

A(J,K).P(J-1,K)  +  B(J,K).P(J.K)  +  C0,K)  PO+TK)  =  G(J,K) 

The  A,  B  and  C  terms  contain  only  DX2,  DX3,  and  p(J,K),  and  the  G  term  contains 
velocities,  viscosities,  and  pressures  for  cells  J,K+1  and  J,K-1.  Similar  equations 
are  obtained  for  all  the  other  cells,  giving  a  set  of  equations  of  a  tridiagonal  form 
which  can  easily  be  solved  by  standard  methods. 


In  the  particular  iterative  solution  method  used  here  the  pressures  are  solved  for 
each  row  of  constant  K  in  turn.  The  velocities,  viscosities,  and  pressures  for 
values  of  K  different  from  the  value  under  consideration  are  taken  from  the 
previous  iteration. 


2.3  Method  of  Solution  of  Equations  for  Velocity  Distribution. 

There  are  2  methods  of  deriving  the  velocity  distribution  based  on  the  model 
proposed  in  Section  2.1.  The  velocity  distribution  within  a  cell  is  parabolic,  as 
given  by  equations  3  and  4.  For  the  sake  of  discussion,  the  velocity  in  the  2 
direction  will  be  considered.  At  the  boundaries  of  the  cells  in 
the  3  direction  the  shear  stresses  in  each  cell  have  to  be  equal.  The  shear  stress  is 
given  by:- 


3i>2 

where  the  lower  case  "v"  denotes  the  continuously  variable  velocity  in  the  cell. 
For  the  boundary  between  cells  (J,K)  and  (j,K-l )  the  equality  of  shear  stresses 
gives:- 


dv2l(l,K) 
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where  v2,G-K)  is  the  value  of  the  velocity  at  the  upper  right  hand  corner  of  the 
(J,K)  cell,  and  v2b(J,K-l)  is  the  velocity  at  the  lower  right  hand  corner  of  cell  (J,K-1), 
see  figure  2.  The  derivatives  of  velocity  are  obtained  by  integrating  equation  3. 

For  all  values  of  K  in  a  vertical  column  of  cells  in  the  channel,  equation  13  leads  to 
a  set  of  simultaneous  equations  in  the  velocities  at  the  cell  corners.  The  values  of 
the  velocities  at  the  channel  walls  are  known  from  the  boundary  conditions,  so  the 
equations  can  be  solved  by  direct  substitution  methods  to  give  the  velocity 
distribution. 

The  same  equations  for  the  velocity  distributions  can  also  be  obtained  in  a  more 
formal  manner  by  converting  equations  1  and  2  into  finite  difference  equations. 

To  calculate  the  velocities  in  the  2  direction,  equation  1  can  be  rearranged  and 
converted  to  finite  differences  to  give:- 


r 


1  V2<J,K+l)-V2(],K) 

ras-(^+ 1)  PX3 

V2(/yK)-V2(/>K-l)N 
-M2UAJ-  DX3  J" 

P(f+1,K) + P(J+l,K+l)-  P(I,K)  -  P(/,K+1) 
DX  2.2 


DX2'(^3^,+^'^' 


V2(I+l,K)-V2(l,K) 

DX2 


MMM 

-n3(/,X).  DX2 


(14) 


where  V2(J,K)  is  the  velocity  at  the  lower  right  hand  corner  of  cell  (J,K).  This 
method  of  differencing  the  equations  was  chosen  because  the  viscosity  has  been 
assumed  constant  in  each  cell,  so  the  local  spatial  derivatives  of  the  viscosity  are 
zero.  Hence  the  derivative  of  the  product  of  the  viscosity  and  the  first  derivative 
of  velocity  can  be  approximated  by  the  difference  of  the  product  between  2  cells, 
and  so  equation  13  should  be  a  good  approximation  to  the  true  equation. 

Multiplying  by  DX3.DX3  and  rearranging  equation  14  gives:- 
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Similar  equations  are  obtained  for  the  other  cells  and  together  they  form  a  system 
of  equations  of  the  form:- 

A(J,K).V2(J,K-1)  +  B(J,K).V2(J,K)  +  C(J,K).V2(J,K+1)  =  G(],K) 

with  the  variables  to  be  solved  V2(J,K+1),  V2(J,K),  V2(J,K-1).  The  equations  can  be 
solved  stripwise  for  constant  values  of  J.  In  this  case  tire  viscosities,  pressures, 
and  V2(J+1)  and  V2(J-1,K)  are  evaluated  from  values  from  the  previous  iteration. 
The  equations  for  each  J  are  then  solved  by  a  direct  substitution  method,  and  the 
process  iterated  until  the  values  of  velocity  converge. 

A  similar  procedure  is  used  to  solve  the  velocities  in  the  3  direction. 


3.  Results  and  Discussion 

A  computer  program  to  implement  the  model  described  above  has  been  written  in 
BASIC.  The  kneading  disc  modelled  had  a  flight  clearance  of  1  mm,  a  depth  in 
the  centre  of  the  channel  of  8  mm  and  a  total  length  of  67.5  mm.  This  channel  was 
represented  on  a  grid  of  cells  numbering  25  by  16  of  dimension  2.7  mm  by 
0.5  mm  (see  Fig.  1).  The  power  law  viscosity  parameters  were  k  =  5000  Pa.sn  and 
n  =  0.5. 

Calculated  values  of  the  pressures  at  the  cell  centres  are  given  in  Table  1.  Only 
the  left  half  of  the  channel  is  shown  because  the  flow  is  symmetrical.  The  positive 
2  direction  is  from  left  to  right  and  the  positive  3  direction  is  down  the  page.  The 
pressure  is  defined  to  be  zero  at  the  centre  of  the  flight  tip  and  at  the  centre  of  the 
channel.  The  values  are  negative  in  the  rest  of  the  region  because  the  flow  is  out 
of  left  half  into  the  right  half  of  the  channel,  where  the  pressure  is  positive.  The 
pressure  increases  with  increasing  distance  in  the  2  and  3  directions  in  the  main 
part  of  the  channel  as  would  be  expected. 

Values  of  the  velocity  in  the  2  and  3  directions  at  the  right  hand  corner  of  the 
cells  are  given  in  Tables  2  and  3.  The  velocities  are  symmetrical  about  the  cells  in 
columns  0  and  12.  The  pressure  gradient  in  the  3  direction  causes  the  upward 
velocity  in  the  3  direction  in  the  left  hand  part  of  the  channel  which  leads  to  a 
circulating  flow.  The  circulation  can  be  seen  in  Tables  4  and  5  which  give  the 
magnitude  and  direction  of  the  total  flow. 

The  FAN  model,  which  was  used  previously  to  calculate  flow  in  a  similar 
channel,  was  essentially  only  a  1  dimensional  model,  and  only  the  variation  of 
pressure  and  velocity  in  the  2  direction  were  calculated  [1].  To  compare  the 
models,  the  FAN  model  was  used  with  the  geometry  and  viscosity  parameters  of 
the  current  model.  The  results  of  the  2  models  are  compared  in  figures  3  and  4. 
Figure  3  shows  the  variation  of  pressure  in  the  2  direction,  and  the  FAN  values 
are  about  10  to  15%  lower  than  the  full  2  dimensional  values.  Given  the 
approximations  involved,  particularly  in  the  FAN  model,  the  agreement  is 
acceptable. 

A  comparison  of  the  calculated  values  of  the  components  of  velocity  in  the  2 
direction  is  given  in  Figure  4.  The  velocities  are  evaluated  along  the  vertical  line 
of  symmetry  in  the  centre  of  the  channel.  It  can  be  seen  that  the  agreement  is  very 
good. 


Pressure(kPa) 


The  model  has  been  constructed  with  the  aim  of  future  expansion  to  3 
dimensions.  The  results  obtained  confirm  the  reliability  of  the  model  in  2 
dimensions  and  do  not  indicate  any  impediment  to  further  expansion. 
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Figure  3:  Comparison  of  calculated  pressure  profile  in  the  2  direction  from  the  current 
2  dimensional  model  and  the  previous  FAN  model. 
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Figure  4:  Comparison  of  the  2  component  of  velocity  in  the  .3  direction  from  the  current 
2  dimensional  model  and  the  previous  FAN  model. 
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Table  1:  Pressure  distribution  (kPa) 


0 

1 

2  3  4  5  6  7  8 

9 

10 

11 

12 

13 

0 

1 

0 

-244 

-535  -514  -429  -345  -275  -219  -173 

-134 

-98 

-64 

-32 

-0 

2 

-0 

-245 

-538  -516  -430  -346  -276  -220  -174 

-135 

-98 

-64 

-32 

-0 

3 

-541  -518  -431  -346  -277  -221  -175 

-135 

-99 

-65 

-32 

-0 

4 

-519  -430  -345  -277  -221  -175 

-135 

-99 

-65 

-32 

-0 

5 

-519  -428  -343  -275  -220  -175 

-135 

-99 

-65 

-32 

-0 

6 

-426  -338  -272  -219  -174 

-135 

-99 

-65 

-32 

-0 

7 

-426  -335  -268  -216  -173 

-134 

-98 

-65 

-32 

-0 

8 

-333  -265  -213  -171 

-133 

-98 

-64 

-32 

-0 

9 

-334  -263  -209  -168 

-131 

-97 

-64 

-31 

-0 

10 

-261  -206  -167 

-129 

-96 

-63 

-31 

-0 

11 

-263  -206  -168 

-129 

-96 

-61 

-30 

-0 

12 

-204  -168 

-129 

-99 

-59 

-30 

-0 

13 

-205  -168 

-128 

-100 

-59 

-31 

-0 

14 

-168 

-128 

-101 

-58 

-31 

-0 

15 

-127 

-101 

-57 

-31 

-0 

16 

-56 

-31 

-0 

Table  2: 

Velocity  distribution  (mm/ s) 

-  2  direction 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

0 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

1 

194 

192 

155 

137 

137 

144 

154 

163 

170 

174 

176 

178 

179 

179 

2 

0 

0 

52 

57 

62 

72 

85 

97 

106 

112 

116 

119 

121 

121 

3 

0 

14 

19 

26 

37 

49 

58 

65 

69 

72 

74 

74 

4 

2 

1 

2 

8 

17 

24 

30 

34 

37 

39 

39 

5 

0 

-3 

-7 

-7 

-3 

2 

6 

9 

11 

13 

13 

6 

-3 

-10 

-13 

-14 

-12 

-10 

-8 

-7 

-6 

-6 

7 

0 

-9 

-15 

-18 

-19 

-19 

-19 

-18 

-17 

-17 

8 

-7 

-15 

-19 

-22 

-23 

-24 

-24 

-24 

-24 

9 

0 

-14 

-19 

-23 

-25 

-27 

-27 

-27 

-27 

10 

-10 

-18 

-23 

-25 

-27 

-27 

-28 

-28 

11 

0 

-14 

-22 

-25 

-27 

-27 

-28 

-28 

12 

-9 

-18 

-23 

-26 

-27 

-28 

-28 

13 

0 

-10 

-18 

-22 

-24 

-26 

-26 

14 

0 

-9 

-14 

-19 

-21 

-21 

15 

0 

0 

-9 

-12 

-12 

16 

0 

0 

0 

4 


17 


Table  3:  Velocity  distribution  (mm/s)  -  3  direction 


0 

1 

2  3 

4  5  6 

7 

8 

9  10 

11 

12 

0 

0.0 

0.0 

0.0 

0.0  0.0 

0.0  0.0  0.0 

0.0 

0.0 

0.0  0.0 

0.0 

0.0 

1 

0.3 

1.1 

2.4 

1.7  -0.1 

-1.1  -1.1  -0.8 

-0.5 

-0.3 

-0.2  -0.2 

-0.1 

0.1 

2 

0.0 

0.0 

2.2 

1.9  -0.9 

-3.0  -3.3  -2.6 

-1.7 

-1.0 

-0.7  -0.5 

-0.2 

0.2 

3 

0.0 

0.6  -2.1 

-5.1  -5.6  -4.6 

-3.2 

-2.0 

-1.4  -1.0 

-0.3 

0.3 

4 

0.0  -2.5 

-6.2  -7.5  -6.5 

-4.6 

-3.0 

-2.1  -1.4 

-0.5 

0.5 

5 

-1.9 

-6.2  -8.3  -7.7 

-5.7 

-3.8 

-2.6  -1.9 

-0.7 

0.7 

6 

-1.4 

-5.4  -8.2  -8..“ 

-6.3 

-4.3 

-3.1  -2.2 

-0.8 

0.8 

7 

0.0 

-4.3  -7.6  -7.9 

-6.4 

-4.5 

-3.3  -2.5 

-0.8 

0.8 

8 

0.0 

-2.9  -6.8  -7.4 

-6.1 

-4.4 

-3.3  -2.6 

-0.9 

0.9 

9 

0.0  -5.3  -6.8 

-5.6 

-4.1 

-3.2  -2.7 

-0.9 

0.9 

10 

-3.2  -5.8 

-5.2 

-3.7 

-3.0  -2.6 

-0.8 

0.8 

11 

0.0  -4.0 

-4.5 

-3.3 

-2.8  -2.5 

-0.8 

0.8 

12 

-2.2 

-3.3 

-2.8 

-2.5  -2.4 

-0.7 

0.7 

13 

0.0 

-1.7 

-2.0 

-2.0  -2.1 

-0.6 

0.6 

14 

0.0 

-0.9 

-1.2  -1.6 

-0.4 

0.4 

15 

0.0 

0.0  -0.7 

-0.2 

0.2 

16 

0.0 

0.0 

0.0 

Table  4:  Magnitude  of  velocity  distribution  (mm/s). 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

0 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

251 

1 

194 

192 

155 

137 

137 

144 

154 

163 

170 

174 

176 

178 

179 

179 

2 

0 

0 

52 

57 

62 

72 

85 

97 

106 

112 

116 

119 

121 

121 

3 

0 

14 

19 

27 

38 

49 

58 

65 

69 

72 

74 

74 

4 

2 

3 

7 

11 

18 

25 

30 

34 

37 

39 

39 

5 

0 

4 

10 

11 

8 

6 

7 

9 

11 

13 

13 

6 

4 

11 

16 

16 

14 

11 

9 

7 

6 

6 

7 

0 

10 

17 

20 

20 

20 

19 

18 

17 

17 

8 

8 

17 

21 

23 

24 

25 

24 

24 

24 

9 

0 

15 

20 

24 

25 

27 

27 

27 

27 

10 

10 

19 

24 

25 

27 

27 

28 

28 

11 

0 

15 

22 

25 

27 

27 

28 

28 

12 

10 

19 

23 

26 

27 

28 

28 

13 

0 

10 

18 

23 

24 

26 

26 

14 

0 

9 

14 

19 

21 

21 

15 

0 

0 

9 

12 

12 

16  0  0  0 
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Table  5:  Direction  of  velocity  distribution  (degrees) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

360 

360 

359 

359 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

358 

358 

1 

2 

2 

2 

1 

1 

0 

0 

0 

0 

3 

0 

358 

6 

11 

9 

5 

3 

2 

1 

1 

0 

360 

4 

1 

67 

69 

42 

21 

11 

6 

3 

2 

1 

359 

5 

0 

151 

140 

130 

112 

74 

34 

17 

10 

3 

357 

€ 

158 

150 

149 

149 

153 

157 

160 

161 

172 

188 

7 

0 

155 

153 

156 

162 

167 

170 

172 

177 

183 

8 

159 

156 

159 

165 

169 

172 

174 

178 

182 

9 

0 

158 

161 

166 

171 

173 

174 

178 

182 

10 

162 

162 

167 

172 

174 

175 

178 

182 

11 

0 

164 

168 

172 

174 

175 

178 

182 

12 

167 

170 

173 

174 

175 

179 

181 

13 

0 

171 

173 

175 

175 

179 

181 

14 

0 

174 

175 

175 

179 

181 

15 

0 

0 

175 

179 

181 

16 

0 

0 

0 

4.  Conclusions 

A  model  of  flow  in  a  2  dimensional  channel  has  been  developed  which  explicitly 
calculates  the  distribution  of  pressure  and  velocity  in  2  dimensions.  The  pressure 
and  velocity  distributions  obtained  appear  to  be  realistic,  and  compare  well  to 
distributions  obtained  with  the  simpler  FAN  method. 

The  model  should  be  capable  of  being  readily  expandable  to  model  3 
dimensional  flows. 
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